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ABSTRACT 


A high-current, high-power X-ray Free Electron Laser (FEL) using the Stanford 
linear accelerator (SLAG) as an electron beam source has been proposed. Such a 
system will provide the scientific community with its first realistic X-ray laser. One 
dramatic use a hard X-ray laser light at the proposed wavelength of 1.5 Angstrom 
would be to image DNA base-pairs [7]. 

With an accelerator already in place, the major cost of developing the SLAG 
system would be in the manufacturing of the undulator. The cost in building an 
undulator is proportional to its length. This thesis uses numerical simulation to 
evaluate the effectiveness of an FEL klystron in order to decrease the length of the 
undulator with a dispersive section. 

Simulations show that the quality of the electron beam injected into the undulator 
can greatly effect the ability of an FEL to produce coherent light, and is a factor in 
determining the length of an undulator needed for saturation. In the high-gain FEL 
(such as proposed at SLAG), a dispersive section is sensitive to energy spread and 
may adversely affect the gain of the system. For moderate to low energy spread, it is 
found that a klystron is useful in shortening the undulator length. 

By examining how the strength of a klystron dispersive section affects the gain of 
a high-current system, we determine a critical strength above which the klystron can 
be detrimental. A dispersive section that is either too strong or not strong enough will 
result in less than optimum gain, or an increase in the required length of an undulator 

Single-mode, phase-space simulations are used to investigate the effect on 
electron bunching and the onset of saturation. Longitudinal multimode simulations 
show the resulting coherence development. 
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I. INTRODUCTION 


With the present world turmoil, the political and economic structure of current 
nations has made it easier for the Third World nations or even small factions of 
terrorist organizations to obtain state-of-the-art weapons. The advances in technology 
as well as hardening of modern day missiles has driven present day ship self-defense 
systems obsolete. The time required to track, engage, and ensure destruction of 
incoming weapons has decreased to seconds. Present day defense systems can no 
longer protect ships adequately. 

The present day ship self-defense system is a layered system that relies on long 
and medium range missiles for long distance engagements. A ship’s gun and short 
range missiles such as RAM and SEA SPARROW are for missiles that get through the 
long range defenses and are at a medium range. The final and most crucial system in 
the layered defense is the Phalanx Close-ln-Weapon System (CIWS), for close in 
missiles threats. All of these systems except for CIWS, are not always available to 
every platform. Even if the outer defense systems are available, the detection and 
engagement of these incoming missiles at long ranges have become difficult with the 
design of sea skimming and low, slow flying missiles. The present self-defense 
systems must also overcome a large scale attack with multiple weapons being fired 
upon them simultaneously. A ship's ability to protect itself against such an attack is at 
present unrealistic. 

CIWS has a reliable detection range of about 5000 m and an effective 
engagement range of less than 2000 m. A typical missile requires numerous hits from 
CIWS before it is destroyed. With the missiles traveling at speeds up to and greater 
than mach 3, the range that CIWS actually kills the missile is much closer than the 
2000 m engagement range. Typically, the actual kill range is around 600 m. If CIWS 
is able to destroy the missile, fragments from the missile are still large enough and 
traveling with enough momentum to reach the ship and create significant damage. If 
the incoming missile is designed to maneuver once it is in close range of its target, a 
CIWS would have to slew its gun barrels causing the rounds to be off target until the 
missile stabilizes and the CIWS is locked on again The time that a CIWS round 
travels to reach the incoming missile also adds to the decreasing range problem. It is 
evident the technology of missiles has outgrown ships’ self-defense systems. 




There is need for a speed-of-light weapon which can reach farther than CIWS 
and kill an incoming missile faster, so as to allow for multiple threat engagements. To 
kill an incoming missile with a laser, the power required on the surface of the missile is 
dictated by the Fluence, F = P\tlA [1]. Fluence (in Jlcm^), is the common measure 
of "hardness" of a missile. It describes the amount of energy absorbed by the 
missile’s skin. A moderately hardened missile might require a fluence of 10 kJIcm^ 
[1]. A spot size A = 5.0 cm^ is achieved at a range ff = 5.0 km, wavelength 
A, = 3.8 microm , with a laser aperture of D = 1.5 meters, where 

._d2j2 

(1-1) 

P is the power received at the missile and Af is the dwell time required to kill the 
missile. For a 0.1 second dwell time, the totetl power to destroy the missile is 
approximately P = 500 kW. The ship’s weapon must emit greater power than this to 
overcome aerosols at sea level over the distance it travels. The power transmitted 
through the atmosphere decreases exponentially as P(z) = Poe"’^, where 
a = 5 X 10"® /cm”’ is the extinction coefficient for a wavelength of X = 3.8 microns, Pq 
is the initial power at the source, and z is the range to the target. At a range of 20 km 
the power required at the source is about 1.4 MW or greater. This amount of power 
has been achieved with chemical lasers. The minimal amount of time, ranging from 
microseconds to seconds, that the missile will travel while the laser is locked on, 
greatly decreases the closing distance the missile travels to only meters. This is a 
significant difference compared to the distance the missile closes in on its target 
before a CIWS obtains a kill, and it is performed at a range that is much greater than 
the engagement range of CIWS. The laser’s speed to target also helps to reduce the 
closure distance. The range at which the missile is hit. even if not destroyed but 
damaged, is sufficient to prevent the incoming weapon from ever reaching it’s target. 

A Free Electron Laser (FEL) has many advantages over conventional lasers in 
producing high power outputs. Conventional lasers of sufficient power use dangerous 
chemicals that are difficult to store and are partially depleted with every firing of the 
laser. The depletion of the source for lasing then limits the number of times you can 
fire, and the number of incoming missiles that can be engaged. An FEL’s power 
source is the electrical power from the ship, which is degraded only if the ship is 
already badly damaged. The FEL is tunable over a wide band of frequencies. That 
allows you to compensate for atmospheric conditions and concentrate the light spot on 
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the inbound missile. 

The Free Electron Laser has other significant uses that can be exploited by the 
military, other than ship’s self-defense. There are many uses in medicine such as 
DNA mapping and X-rays. Target tracking and identifying can be enhanced with a 
tunable laser. Ideally, the burden of research and development can be shared by 
civilian and military ensuring continuous upgrading and off-the-shelf products. 
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II. BACKGROUND 


A. FEL THEORY 

There are two major requirements for a FEL to produce emission; a relativistic 
electron beam supplied by an accelerator, and an undulator with a spatial periodic 
magnetic field [2]. The FEL creates coherent, monochromatic light by accelerating 
relativistic electrons in the transverse direction. The periodic magnetic field from the 
undulator exerts a force on the relativistic electron beam as it passes through the 
undulator causing the electrons to "wiggle" in the transverse direction. The electrons 
are accelerated from side to side causing them to spontaneously emit light in a 
forward cone along the electron beam axis. To collect a fraction of the spontaneously 
radiated electromagnetic field, mirrors are placed on the beam axis beyond both ends 
of the undulator creating an optical cavity. The coupling created by passing the 
electron beam through the stored electromagnetic field causes stimulated emission. 
As a result, the optical field grows and coherent light is produced from the FEL [2]. 

The electrons are traveling relativisfcally with energy ymc^ prior to entering the 
undulator, where m is the electron rest mass and c is the velocity of light. The 
Lorentz factor of the electron is y= 1 W 1 -P®, where p = v/c and v is the magnitude of 
the electron velocity. Upon entering the undulator, a force from the magnetic field acts 
upon the electrons causing the electrons to move back and forth in the transverse 
direction. 

The optical field properties can be controlled by controlling the parameters of the 
undulator An undulator can be of different designs using either permanent magnets 
or electromagnets with either linearly or circularly polarized magnetic fields producing 
linearly or circularly polarized light. The strength of the magnetic field determines the 
magnitude of deflection of the electrons and therefore the electron z velocity. Varying 
the electron z velocity (i.e. varying the electron energy) makes an FEL tunable over a 
large range of wavelengths. The undulator parameter, K = eBx^Znmc^, where B is 
the root-mean-square (rms) value of the undulator magnetic field, characterizes the 
undulator’s magnetic field strength. In a linear undulator B = S/Vl, because of non¬ 
constant acceleration, while in a helical undulator B = B, because of the continuous 
acceleration, where B is the peak magnitude of the undulator magnetic field. 
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While the quantum mechanical approach shows that the emission of a photon 
causes a recoil in the electron, a classical approach to the FEL theory neglecting recoil 
can be used because the energy loss is small. The ratio of energy lost with electron 
initial energy is on the order of 10"''. It requires many such emissions to create an 
effect that would be significant, due to the high energy of the electrons. With the 
optical beam, the photon density is so great that the statistical fluctuations of the 
optical beam can also be ignored [11]. 

The classical approach describes an electron wiggling and emitting light while 
transversing the undulator’s magnetic fields. The light is emitted along the electron 
beam axis into a fonward cone. Figure 2.1 [4] shows the interaction between electrons 
and the optical field. It shows that the optical electric field acting on an electron will 
take energy away from an electron. To accomplish this transfer of energy the optical 
field and the electron must have the correct relative phase. The spread of electron 
phases in comparison to the optical field causes some of the electrons to take energy 
from the optical field when they are out of phase. The electrons then form bunches 
over the optical wavelength. The end result is that coherent light is formed from these 
electron bunches. 


As mentioned, unlike classical lasers where coherent light is produced by 
electrons transitioning from one quantum state to another, the FEL produces light by 
the interaction of a relativistic electron beam and external forces. The FEL is then 
able to generate coherent light at tunable frequencies depending on the undulator 
parameter K. As the electrons are relativistic the, undulator wavelength Xq is Lorentz 
contracted to A,'o = in the electron beam frame. The optical field passing over the 
electrons is Doppler shifted to a wavelength X' = (1 + = 2yX, where X is the 

optical field wavelength in the lab frame [3]. To have sufficient energy exchange 
between the electron and the optical field a condition of resonance must occur, 
X' = Xq' in the beam frame. In the undulator frame 


X = 


Xo(1 + K^) 

2f 


(2.1) 


By evaluating Eq. 2.1, it is apparent that an FEL can be continuously tuned by 
changing either the electron energy ymc^, the undulator’s magnetic field strength B, or 
the undulators wavelength Xq. The undulator can also be designed to compensate for 
the falling out of resonance by the electrons while they give up energy. By changing 
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Figure 2-1. Fundamental physics of the FEL interaction. 



















>^0 through the undulator the electron can be kept in phase with the optical tield longer. 
This is called tapering the undulator, which increases the efficiency of the FEL. 


B. ELECTRON DYNAMICS AND THE PENDULUM EQUATION 

The individual electron’s motion inside an undulator is described by solving the 
classical equations of motion. The Lorentz force is the only significant force acting on 
an electron traveling down the length of the undulator [3]. In the presence of the 
undulator and optical fields, the Lorentz force equations are 



(2.1a) 


(2.1b) 

Y-2 = 1 - (^ • ^ . 

(2.1c) 

where ^ is the electron normalized velocity vector, 

helical undulator is circularly polarized and is given by 

The magnetic field of a 

Su = e [cos{kQz), sinikoz), 0] , 

(2.2) 


where the undulator wave number = 2nlXo. and z is the direction of the undulator 
axis. The optical electnc and magnetic fields are assumed to be circularly-polarized 
plane waves and are given by 

= E [cos{\v),-sin{\\i), 0] , = E [sin{\\i), cos{^i), 0] , (2.3) 

where E is the optical field magnitude and \); = kz + wf -k)). The optical wave number 
is k = 2nlX, the optical frequency is <o = ck, and is the phase of the optical wave. 

Substituting the undulator magnetic field Eq. 2.2 and the optical wave 
electromagnetic fields Eq. 2.3 into the Lorentz force equations Eq. 2.1 , they become: 

“ [E(2-p^)(cos(v|;).-s/r}(vi;),0) + B[-sin{koz). cos(koz). 0)](2.4a) 
^z) = ^ [£(PxCos(m/) - PySin(w)) + B (PxS/n(koz), - p^cos(V)] (2.4b) 
^=--^E [P;,cos(m/) - p^sin(M/)] , (2.4c) 

where = (P;,, Py, 0). As the electrons travel down the undulator in the z direction 



they are acted upon by a transverse force from the optical electric and magnetic fields 
and the magnetic field from the undulator. The force from the optical field is 
proportional to E (1 - while the force from the undulator magnetic field is 

proportional to B. For a relativistic electron p^ 1, so 1 - p^ -> 0. Therefore, we 
can assume E (1 - p^) c S p^. The force acting on the electron in the transverse 
direction is mainly due to the undulator magnetic field. Equation 2.4a can then be 
approximated as 


Pz B[-sin(koZ), cos{koz). 0] . 


(2.5) 


On integrating Eq. 2.5, the electron’s velocity in the transverse direction is found to be 


^4, = -^[cos{koz). sin(/coz). 0] 


(2.6) 


where K = eBkf,l2nmc^ = eBlkomc^ is the undulator parameter. 

Examining the transverse and longitudinal motion of the electron as it travels 
down the undulator, it is noted that the transverse motion is small, p^ = K/y >*: 1. 
Inserting Eq. 2.7 into Eq. 2.5c gives the electron energy evolution. 



(2.7) 


where the dimensionless electron phase in the combined undulator and optical fields is 
^ + (j) = + koz = {k + k^z - cot + (t>. The initial value of the electron phase is 
C(0) = Co = (^ + ^o)^o- Since k fhen Co “ = SrtZo/A.. when y »> 1. The 

electron phase C determines the electron’s z position relative to a single optical 
wavelength X. By using Eq. 2.1c and inserting in Eq. 2.6 gives 


y-* = 1-P,2-P,^ = 1-P,2-4 . 

r 


(2.8) 


(1 + K^Yr^ = 1 - p,2 . 


(2.9) 


Then, take the time derivative of both sides 


-2y-3y (1+K^ = -2p,p^ 


(2.10) 


which relates y to p^ as 


y_ 

Y (1 + K^) ■ 
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(2.11) 


(2.12) 


Taking the first derivative of the electron phase ? gives 

t = (k + ko)z -<o = {k + ko)c^^ - co . 
and a second derivative gives 

t = {k+ko)c^, . 


By using Eq. 2.11, and Eq. 2.13, and the resonance condition 

X =_ fi « A_ 

Y (1 + K^)(k + ko) <0(1 + K^) 2o)o ■ 
where wq = co{1 + K^)l2:f. The electron’s equation of motion Eq. 2.7 
with Eq. 2.14 to give 


c 


2(OoeKE 

^mc 


cos(C + 


(2.13) 


(2.14) 
be combined 


(2.15) 


Define the dimensionless time x = ctIL where L is the length of the undulator to 
make the equation of motion, Eq. 2.15, dimensionless. As an electron passes through 
the undulator. x = 0-> 1. The dimensionless optical field is defined as a = lale'*, 
where la I = AnNeKLEly^mc^, and N is the number of periods in the undulator. The 
pendulum equation can be written as 

F = v = lal cos(C + <ti) , (2.16) 

where C is the second derivative with respect to x. v is the first derivative with respect 
to X, and the dimensionless electron phase velocity is defined as 
V = L [(k + /ro)Pz - k]. When v = 0, the electrons are at resonance with the optical 
wave and X = Xo(1 + Changes in the electrons’ energy correspond to giving 

or taking energy from the optical wave. As the electron decreases its energy, its 
phase velocity decreases and it gives up energy to the optical wave. The maximum 
energy loss of the electron occurs when ^ + <l> = n- 


C. THE FREE ELECTRON LASER WAVE EQUATION 


The development of a complex optical wave equation starts from Maxwell’s wave 
equation acting on a vector potential X(f.f) driven by the current density 





(2.17) 
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where the vector potential X(z.f) is 

X{z.t) = E(z.t) [sin(;cz - cof + (t)(z,f)). cos(kz - wf + (t)(z,0), 0] , (2.18) 

and the electric field ^ is 

f = E(z,f)[cos(kz - cof + <t)(z.f)),-sin(fcz - cof + 4i(z,f)), 0] , (2.19) 




C dt 


(2.20) 


Since the wave is traveling only in the z direction the Maxwell wave equation 
operator then is 

__1 ii 

dz^ af2 ■ 

Taking the first derivative of the vector potential X with respect to z gives 


( 2 . 21 ) 


[sin(M/), cos(h/), 0] + -f |k + || [cos(m/), -sin(M/), 


dz k dz 


0] , 


(2.22) 


where vy = fez - cof + (!>, and c/co = 1/#e. Taking the second derivative with respect to z 


0-if 

+ J [cos(m/), -sin(v|/). 0] + j(k + ||)Vsin(M/), -cosM, 0] . (2.23) 


Assuming slowly varying amplitudes and phases where dE/dz « kE, d^ldz <*: k^, 
dEldt < <o£, d^ldt < co<^, and (o* kc allows us to use only first-order terms and 
simplify Eq. (2.23) to 

0 = 2 If [cos(m/). -sin(v). 0] + f ['f + H] t-sin(v). -<os{w). 0] . (2.24) 


+ =fe- + 


(2.25) 



knowing {d^ldzf -« d^ldz, and putting Eq. 2.25 into Eq. 2.24, it becomes 

^ = 2 [cos(v)/). -sin(M/). 0]+k£+2E|^ l-sinM. -cos(h/). 0] . (2.26) 

dz‘^ oz [ dzj 

The time derivative portion of the Maxwell’s wave equation operator is 

^ “s(h'). 0] . (2.27) 

Taking the first derivative with respect to time gives 

1^ [E (sin(H/), cos(n/), 0)] = ^[sin(M/). cos(w). 0] 

+ E - <o| [cos(vv). -sin(v,/). 0] . (2.28) 

The second derivative gives 

^[sin(vt/), cos(m/). 0] = - (ojcos(M/). -sin(H,). 0] + 0[sin(vt/), cos(vi,), 0] 

^ ^[if " £||[cos(m/). ^in(vi,), 0] 

+ E - (oj [-sin(H/). ^os(M/), 0] . (2.29) 

Again, assuming the slowly varying amplitudes and phases allows us to neglect all but 
first-order terms, so Eq. 2.29 becomes 

^ = -2 c ^[cos(m/). -sin(M/), 0] -r E - o)j [-sin(M/), -cos(vy), 0] (2.30) 
By expanding 

again knowing (3<t»/3f)^ ■« d^ldt, we put Eq. 2.31 into Eq. 2.30 and it becomes 

^ = -2 c ^[cos(h,), -sin(M/), 0] -r c (<o E - 2 E |^)[-sin(v,/), -cos(m/), 0] (2.32) 
Putting Eq. 2.26 and Eq. 2.32 into Eq. 2.17 gives 
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_ JL 

|X = 2[^.1 



1 .^1 

[dz^ dt^j 

[3z c 

afj ’ 

[az 

c at J 


(2.33a) 


where the unit vectors are = [cos(\)/), sin(\)/), 0] and = [-sin(4;), -cos(v)/), 0]. We 
also know Eq. 2.33a is a multiple of the current density in Eq. 2.17 so 


471 

c 




2 


1 

c dt\ 


. 


(2.33b) 


To solve for the time variation of the field envelope make the change of variables 
z ^ 2 - c t. Then 2 = z + c t where the undulator length L = NXq, N is the number 
of periods, Xq is the wavelength of the undulator, and x = ctiL is the dimensionless 
time. By making the change of variables, we change the partial derivatives. 


_^ ^ ^ ^_.j 

dz ~ dz d2 d2 dr ~ B2 


(2.34) 


and 


^_^ ^ ^ ^ _^ ^ J_^ 

C dt C dt d2 C dt dr d2 L dx ' 

Putting Eq. 2.34 and Eq. 2.35 together simplifies the differential operator to 
dz C dt L dx ' 


(2.35) 


(2.36) 


C ^ L dx 

Doing the same to Eq. 2.33b but projecting onto ^ gives 


(2.38) 


The current density for a single electron is = - ec^^. The electron trajectory 
will be found to be - {Kly) [cos(koz), sin(/Coz), 0], so that 

[cos(koz), sin(koz), 0] , (2.39) 

where y= lWl - v^/c^. v is the velocity of the electron, is the undulator wave 
number 2tiIXo, and the undulator parameter is K = eBlmk^c . 
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Inserting Eq. 2.39 and the value for?, into Eq. 2.37 gives 

^ - 25^ [cos(koz). s\n(koz). 0] ■ [cos{n/). sinM, 0] . (2.40) 

or Y 

On performing the dot product, Eq. 2.40 becomes 

||= - 25^cos(C + .t.) , (2.41) 

where the electron phase is C, = (k + /fo)z -<of. 

The wave equation Eq. 2.41 is for a single electron. The average value over 
many electrons gives 



where p is the electron particle density and 

<cos(C + ■!.)>= L 


<cos(C+ <(')> . 

cos(C, + (!>,) 

Ne 


(2.42) 


(2.43) 


Similar manipulations of Eq. 2.38 gives us 

<sin(C H 

ax Y 

The complex electric field is E = IFI e' ♦ and 


aiEi 

ax 


When Eq. 2.42 and Eq. 2.44 are inserted into Eq. 2.45, we get 

IE ^ _ 2nLeKp ^ ^ ^ 


From the pendulum equation, we know that our dimensionless optical wave field 
is a = I a I e' where I a I = AnNKLe IE l/Y®mc®. so that 

^ ^ ° _8jle^pNf^ , , * ,2.48) 


■^mc^ 
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The dimensionless current is j = STj^e^pNK^L^I^mc^. Which leads us to the 
Free Electron Laser optical wave equation 

a = -j <©-' ?> . (2.49) 

Looking at the real and imaginary parts of Eq. 2.49, you have 

la I = -;■ <cos(C + (!>)>, I =-r^ <sin(^ + ((.)> , (2.50) 

lal 

where lal is the optical field amplitude and ^ is the optical field phase. Equation Eq. 
2.50 shows that bunching of electrons in phase at (C + <t>) = n drives the optical wave 
amplitude and produces a gain, while bunching electrons around (C + (|)) = n/2 will 
drive the optical phase. Increasing the electron beam current density ; will increase 
gain and optical phase evolution. The optical fields increase in magnitude will reduce 
the effects of the electron phase evolution. 


D. LOW GAIN FEL 

The Maxwell-Lorentz (standard classical, i.e. non-quantum dynamics) theory for 
an FEL can be formed by combining the simple pendulum equation, Eq. 2.16, and the 
optical wave equation, Eq. 2.49. The strength of the dimensionless current density j 
governs the coupling between electrons and light. It determines the response of the 
optical wave to the bunching of electrons in phase C The Maxwell-Lorentz theory is 
generally valid for both weak (la I « n) or strong (lal » n) optical fields [ 5 ]. 

From the pendulum equation, Eq. 2.16, it is evident that the electron energy loss 
is greatest when the electron bunch is at C + fti = rr. which provides for the maximum 
gain in the optical field. As the electron bunch moves toward C + <!> = 0, over-bunching 
occurs and energy is absorbed back into the electron beam. 

Maximum coupling between the optical field and the electron beam occurs when 
the initial electron phase velocity is at resonance, vq = 0 [6]. A "realistic" beam of 
monoenergetic electrons will enter the undulator with a random spread in phase C. 
There will be just as many electrons gaining energy from, as losing energy to, the 
optical field, resulting in zero gain. To achieve a useful gain, the electrons must have 
a velocity slightly above resonance so that more electrons give up energy than absorb 
it in an exchange with the optical field. 
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To understand the operation of an FEL, it is useful to examine the evolution of 
the electrons as they move through phase space (C,v). The path in phase space for 
each electron is dependent on the initial values of Co vq. A loi«-gain FEL 
produces an electron phase-space path that can be represented by a simple 
mechanical pendulum [3], and is shown analytically by 

v2 = vo^ - 2 la I [sin{C -t- <!)) - sin(Co)] ■ (2.51) 

By examining Eq. 2.51, we can note that there are fixed points, {ti/ 2, 0) and (37i/2, 0), 
where the pendulum does not evolve. The electrons near these phase space points 
evolve slowly. The electrons that are at the stable fixed point (n/2. 0) correspond to a 
pendulum at the bottom of its arc. The electrons with a large phase velocity 
Ivl > 21aare in an open orbit. This would correspond to a pendulum swinging in 
one direction continuously about its pivot point. An electron in a closed orbit is like a 
pendulum swinging with a small amplitude oscillating about the stable fixed point. A 
change in ^ would cause a change in the horizontal position of the pendulums fulcrum. 
The phase-space path that separates the open and closed orbits of an electron is 
called the "separatrix" and is given by 

Vs2 = 2lal[1-sin(Cs-Kj))] . (2.52) 

The separatrix has a peak-to-peak height of 4lal’'^. Its horizontal position depends on 
the optical phase ()). As the optical field grows, so does the separatrix height. 

Since the electron’s evolution in phase space is periodic over each optical 
wavelength X, the interaction in an FEL can be modeled by tracking sample electrons 
over one wavelength. Figure 2-2 is an example of a phase space plot with low current 
y = 1, a moderate optical field strength of la I =4, and an electron phase velocity of 

Vo = 2.6. 

The electrons enter the undulator at t = 0, in a uniform distribution of initial 
phases. Co = - -> 3n/2. The electrons then interact with the optical field of 

strength ag = 4. The initial electrons are displayed as light gray dots. As they evolve 
through the undulator they become darker until they become black dots at the end of 
the undulator. x = 1. The electrons to the left experience an increase in phase velocity 
which causes them to move to the right, while the electrons to the right experience a 
loss of energy and decrease in phase velocity. This produces a spatial bunching 
observed near phase C = n at x = 1. so more electrons are giving up energy than 
taking energy away from the optical field. There is a limit in the process of energy 
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exchange where the electrons have given up so much energy that their phase and 
phase velocity cause the electrons to start absorbing energy back from the optical 
field. An undulator’s length can be designed to prevent the electrons from reaching 
this stage. 

The optical gain is defined as the fractional increase in optical field power over 
time T, 


G(t) 


(a(T)2- ao®) 
a 2 


(2.53) 


Gain is pictorially displayed in the upper right graph of the phase space Figure 2-2. 
The electrons in the beam are assumed to be uniformly distributed in phase Co and 
have the same initial phase velocity vq. A change in electron energy corresponds to a 
change in electron phase velocity. The average energy lost by an electron is 
y/77c^(<v> - vo)/4nA/. The energy lost by the electron determines the optical gain to be 
G = 2 y(vo - <v>)/ao2 [ 3 ], The phase velocity average is found by expanding the 
pendulum equation. Eq. 2.16. The small signal gain equation becomes [3] 


G(t) = j 


\ 2 - 2cos(vo-c) - vqt sin(voT)] 


(2.54) 
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If the phase velocity is far away from resonance, IvqI »n, the final gain at t - 1 will be 
small. Optimum coupling, resulting in maximum gain, occurs when IvqI < n The 
natural gain bandwidth for an FEL is derived from the relation Avq = 4nNAyly. The 
natural gain bandwidth is 


M.\ 

y \ 


2N ■ 


(2.55) 


An FEL gain spectrum and optical phase shift is plotted in Figure 2-3, for a weak 
field ao = 1, low current 7 = 1, and the initial phase velocity between the interval 
-12^vo^12. The gain curve is antisymmetric about vq with maximum gain at 
Vo = 2.6. 



Figure 2-3. Weak-field, low-gain FEL gain spectrum and optical phase shift plot. 





E. ENERGY SPREAD 


The previous section dealt with a perfect beam where the electrons were 
uniformly distributed over the optical wavelength with the same initial energy. This is 
not achievable, since there is a trade-off between electron beam quality and current 
density. Increasing current / tends to decrease beam quality from the accelerator, 
while increasing the length of the undulator by increasing N increases the PEL’S 
sensitivity to beam quality because the natural gain bandwidth is narrower. Where the 
natural gain bandwidth is defined as the range of phase velocities where the growth 
rate is reduced by 10 % below its peak value [3]. 

There are two qualities that describe electron beam quality: the electron beam 
energy spread, Ayly, and emittance, e. Electrons with same initial phase but 
different phase velocities will drift apart on the z axis [6]. At the end of the undulator 
their final phase difference is roughly AC » Avx, where x = 1, for each electron. It is 
apparent, for a change in initial energy Aymc^, there is a corresponding change in 
electron phase velocity, Av « ^nNAyly. A Gaussian distribution of energies with an 
rms spread , Ayly, is characterized by standard deviation oq = AnNAyly. When 
oq = K, bunching becomes impaired [6]. Typical energy spreads from an accelerator 
are on the order of Ayly ~ 0.001. 

Emittance, e, is the average angular spread times the average beam radius, 
E = TQ, where r is the average beam radius and 0 is the average angular divergence 
of the electron beam. The normalized emittance, e„ =ey= r0y, is a measure for beam 
quality. As the electrons are accelerated and beam energy increases, the normalized 
emittance remains constant. There is a Gaussian spread in angles with rms value 0 
resulting in an exponential distribution of phase velocities. The standard deviation in 6 
is oe = 47x^7^ 9^/(1 + K^). As with oq. when cg is equal to or greater than n, electron 
bunching is impaired. 

F. HIGH-CURRENT HIGH-GAIN FEL 

In high current, high gain FELs j n, the magnitude of the optical field changes 
significantly over a single pass of the electrons through the undulator. This is evident 
in the wave equation, Eq. 2.48. As we did in the low gain FEL, assume a perfect 
quality beam operating on resonance, vg = 0. From the pendulum and wave 
equations, the optical field grows exponentially for large j as it travels down the 
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undulator [3]: 

la(T)l = ^ e0/2)''^./2 . (2.56) 

and the phase evolves as; 

f ■ (2.57) 

The PEL gain for large j is 

G(t) « . (2.58) 

For an electron beam that is oft resonance, the single pass gain is 

= , (2.59) 

where xg = {21 is the bunching time. The peak gain will occur slightly above 
resonance. The phase velocity has now become a function of current. As y -> « the 
gain spectrum becomes symmetric centered at resonance. 

The gain bandwidth is Avq = 2^^. The bandwidth for a low gain PEL is 
Avq = 71/2. It is evident that the bandwidth for the high gain case can be much larger 
than that for the low gain PEL. 

G. THE KLYSTRON UNDULATOR 

The purpose of a "dispersive section" in a klystron undulator is to increase 
optical gain in a weak optical field. There are three major sections to a klystron 
undulator. The undulator that was described earlier is divided into two sections 
separated by a ""drift section"" or a "dispersive section". The first section of the 
undulator is called the '"modulator"". It develops the phase velocity differences to 
prepare the electrons for bunching. The next portion is either a drift section, i.e. an 
open section that allows the electrons to continue their bunching without interaction 
with an optical field, or a dispersive section, made up of magnets that force the 
electrons far off resonance. The final section of the undulator is called the "radiator", 
where the electrons and the optical field interact and create the coherent optical beam. 
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Mathematically, there is no difference between the use of a drift section or a 
dispersive section [3], The drift section requires a large distance over which the 
electrons travel with no optical interaction. This is easily achieved, but can possibly 
double the length of the undulator, or even larger. A dispersive section forces the 
electrons off centerline in an arc, and allows for the use of a shorter undulator. The 
strength of either type of sysfem is given by the dimensionless drift time D, where D 
can be determined by the number of undulator periods equivalent to the drift space, 

D = ^. (2.60) 

where /V^ is the length of the drift space in undulator periods. 

While the electrons are traveling in the modulator and radiator the equations 
governing their motion are still Eq. 2.16 and 2.49. While in the dispersive or drift 
section there is no electron-optical wave interaction. The only change in electron 
phase is caused by the constant phase velocity, v; therefore, the change in phase 
velocity is Av = 0, so AC = vD is the only change in electron phase. The electrons 
then will bunch without any interaction from the optical field, and the optical field 
strength will not grow. 

When using a klystron the gain and phase equation are [3] 

G(vo) « ^ sin(voD) . (2.61) 

and 

't>(vo) « ^ cos(voD) . (2.62) 

It is apparent that the maximum gain is about jDi4 when vq = n/2D. It is also 
apparent that with low j, the addition of a klystron enhances gain. The spread in initial 
phase velocity is critical. The natural gain bandwidth becomes narrower, dictating that 
beam quality improve. The spread in initial phase velocity should be kept to 
Avo S nID, which is smaller than that required for a regular undulator. The stronger 
the dispersive section, the greater the need for increased beam quality. 

Hgures 2-4a and 2-4b are the optical gain spectra for a low-current (; =1), 
weak-field (ao = 1), PEL without a klystron and with a klystron, D = 2. Along with the 
narrowing of the natural bandwidth, gain is significantly improved by addition of a drift 
section or dispersive section. 
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Figures 2-5a and 2-5b are the optical gain spectra for a high-current (y = 100), 
weak-field (ao = 1) FEL without a klystron and with a klystron of magnitude D = 2. 
The natural bandwidth is increased and gain is significantly improved with the use of a 
klystron undulator even in a high-current FEL. The size of D for a high-current FEL is 
significantly smaller than that for a low-current FEL. 



Figure 2-4a. Low-current, weak-field FEL gain spectrum, without dispersive section. 
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**** Gain Curve **** 
j=100 aQ=l 

75.67 

0.0 

- 75.67 

Figure 2-5a. High-current, weak-field FEL gain spectrum, without dispersive section. 




Figure 2-5b. High-current, weak-field FEL gain spectrum, with dispersive section. 


23 












24 



III. SLAC X RAY FREE ELECTRON LASER 


A. INTRODUCTION 

Research in the area of high power FEL’s is important to the military for the 
development of missile defense. The research and development that is being 
performed at SLAC is in the high peak power (GeV) FEL regime. By developing hard 
X ray laser light at a wavelength of 1.5 A, imaging of DNA base-pairs may be 
accomplished [7]. The FEL is the best possible source for an X ray laser. 

The SLAC project proposes to develop coherent hard X ray radiation at 
wavelengths as short as 1.5 A. It will use a single pass of an electron beam through 
an undulator with no initial optical field. The optical field will be started from noise 
using a process called Self-Amplified Spontaneous Emission (SASE). There are no 
mirrors for feedback of the optical beam. The electron beam makes a single pass 
through the undulator creating a cone of light in the direction of the electron beam’s 
path. The only interaction between the electron beam and the optical field happens 
during this single pass. Figure 3-1 is a diagram of the system proposed by SLAC. It 
combines the use of the SLAC linear accelerator (linac) with an undulator to produce 
the coherent X ray light. The mirror merely directs the light beam at grazing angles. 



The system design would incorporate a section of the SLAC linac accelerating 
the electron beam to an energy of ymc^ = 15 GeV with peak current of / = 5000 A. 
Table 3-1 contains the design parameters for two systems that would produce either a 
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4.5 A or 1.5 A laser beam [8]. 


PARAMETERS I 

A 

B 

A. (Radiation Wavelength) 

1.5 A 

4.5 A 

Aq (Undulator period) 

2.67 cm 

4.0 cm 

B (Peak undulator field) 

1.173 T 

1.66 T 

Undulator type 

Helical 

Linear 

K Undulator parameter (peak) 

3.4 

4.2 

ymc^ (Electron energy) 

15GeV 

15GeV 

/ (Peak current) 

5000 A 

5000 A 

Ayly (Energy spread) 

0.0002 

0.0002 

L (Undulator length) 

25 m 

38 m 

Output power 

11 GW 

66 GW 


Table 3-1. Proposed parameters for an X ray PEL at SLAG. 

The dimensionless current for / = 5000 A is j = 2056 for the 4.5 A laser and 
y=3121 for the 1.5 A laser. As mentioned earlier, the dimensionless current 
determines the response of the optical wave to the electron beam bunching. The gain 
of the system, because of the high j, is in the high-gain regime. Figure 3-2 is a 
phase-space plot for the 4.5 A system and figure 3-3 is for the 1.5 A system. Each 
simulation is started from an optical field ao = 0 and the electrons are all started with 
random initial phase C between -7t/2 and 3n/2. The simulations ran with an undulator 
length of /. = 35 m. A filling factor F = 0.1 was used, where the filling factor is the 
ratio of the area of the electron beam to the optical beam. The electron beam radius 
r^, =0.0015 cm. The initial electron phase velocity has a standard deviation of 
oq = 2.2 and 4.2, where oq = 4nNAyly. The average initial phase velocity is centered 
around vq = 0. The electrons, in both scenarios, have bunched and are on the 
bottom-half in the separatrix at the end of the undulator, x = 1. The electrons have 
given energy to the optical field. The phase-space simulations indicate the production 
of an X-ray beam with saturation of the system at undulator lengths of 34 and 35 
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meters, respectively. The logarithmic power evolution shown in the upper right-hand 
corners of each figure indicate the systems have reached saturation. The peak power 
achieved for the 4.5 A is approximately 66 GW and for the 1.5 A it is approximately 11 
GW. 



Figure 3-2. Phase-space plot for 4.5 A system. 



Figure 3-3. Phase-space plot for 1.5 A system. 
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In either scenario, the proposed parameters would create a hard X ray FEL as 
indicated by the phase-space simulations. A large portion of the FEL cost is in the 
production of the undulator; reducing the length of the undulator reduces the FEL cost. 
A method of reducing the length of the undulator is by introducing a dispersive section 
in the undulator to make an FEL klystron. A penalty resulting from the klystron design 
is the onset of saturation at weaker fields than without the klystron [9]. 


B. KLYSTRONS AND HIGH-CURRENT HIGH-POWER FEL 


An FEL with an optical wavelength of A. = 40 A was proposed using the SLAC 
linac as an electron beam source [10]. Table 3-2 gives the original proposed 
parameters for the optical wavelength of A. = 40 A. The system required an undulator 
with a length of 60 meters. The phase-space plot in Fig. 3-4 shows a simulation of a 
soft X-ray FEL at 40 A using the parameters in Table 3-2. an initial radiation field 
ao = 0, and a phase velocity spread of standard deviation oq = 3.6 centered around 
Vo = 0. 


PARAMETERS 

VALUE 

A (Radiation Wavelength) 

40 A 

Ao (Undulator period) 

8.3 cm 

K (Undulator parameter) 

4.2 

Undulator type 

Linear 

(Normalized rms emittance) 

3.0 mm-mrad 

ymc^ (Electron energy) 

7.0 GeV 

/ (Peak current) 

2500 A 

Energy spread 

0.0004 

L (Undulator length) 

60 m 

Output power 

28 GW 


Table 3-2. Proposed parameters for a A = 40 A X ray FEL at SLAC. 
(Without dispersive section) 




*** FEL Phase Space Evolution *** 
j=4000 a^=0 N=732 



-ti/2 K 371/2 0 T 

Figure 3-4. Phase-space plot for 40 k system. 
(Without dispersive section) 


The optical field grows from noise as a result of interaction with 200,000 
electrons randomly positioned in the initial electron phase. The dimensionless current 
j = 4000 corresponds to a cument I = 2500 A. Notice, at the end of the undulator, 
t = 1, the electrons have already reached saturation as shown in the power plot in the 
upper right diagram. Having gone through one complete transversal inside the 
separatrix, the bunched electrons are near the bottom of the separatrix, for a second 
time, where maximum energy transfer to the optical beam occurs. When the electrons 
are bunched at the upper half of the separatrix the electron beam has taken back 
energy from the optical beam as indicated by the drop of the optical power after the 
onset of saturation. Reducing the undulator length to approximately 40 m, at the onset 
of saturation, will result in peak power at the end of the undulator. The majority of 
electrons will be at the bottom of the separatrix and the electrons will have given up 
the maximum energy possible to the optical beam. By reducing the length of the 
undulator we reduce the dimensionless current, j = Bn^e^pNK^L^FI'^mc^. With an 
undulator length of 40 m, the dimensionless current is / = 2000. Figure 3-5 is the 
phase-space plot for such a system. 

To reduce costs of manufacturing an undulator of such a long length and 
maintaining the same results, a dispersive section can be added to an undulator, thus 
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making a klystron. Table 3-3 lists the proposals for a system that uses a klystron to 
produce a 40 A PEL [11]. 



The PEL design at X. = 40 A with a klystron has an undulator length of L = 17 
meters, electron beam energy of 5 GeV, and a peak current of 5000 A, which 
corresponds to a dimensionless current density of j = 932. This PEL was examined 
without a klystron and then with a klystron. 

The case without the klystron is shown in Pig. 3-6. The electron bunch is at a 
phase where it gives energy to the optical field at the end of the undulator. The 
dimensionless power P corresponds to a peak power of 15 GW. In Fig. 3-7, the use 
of the klystron PEL with the same parameters and D = 0.5 shows that electrons have 
bunched in the lower half of the separatrlx and at the optimal phase relation ^ + i|i = n 
for maximum gain. The peak power calculated for the klystron PEL is 31 GW. While 
not shown, saturation will occur earlier with the use of a dispersive section than 
without, reducing the overall length of the undulator. 
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PARAMETERS 

VALUE 

X (Fund. Rad. Wavelength) 

40 A 

Xq (Undulator period) 

4.0 cm 

K (Undulator parameter) 

4.2 

Undulator type 

Linear 

En (Normalized rms emittance) 

3.0 mm-mrad 

ymc® (Electron energy) 

5.0 GeV 

/ (Peak cun^ent) 

5000 A 

Energy spread 

0.0004 

L Undulator length 

17 m 

Output power 

2 GW 


Table 3-3. Proposed parameters for a X. = 40 A X ray PEL, 
(With dispersive section). 
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To generalize the optical field in the longitudinal direction, several sites along the 
complex wave envelope. a(z), are followed. The length in the z direction is scaled to 
the slippage distance NX, the distance the electron pulse slips back relative to the 
optical pulse over one pass through the undulator. This is the distance over which the 
optical field and the electron beam interact. All longitudinal distances are normalized 
to this slippage distance so that zINX -» z. This then gives 



Figure 3-7. Phase-space plot for 40 A X ray FEL, 
(With dispersive section). 


a a{z), V v{z), and ^ ^(z) over a number of sites. The pendulum and wave 

equation can be generalized to [3] 

^^ = la„cos(C,„ + 4.,) (3.1) 

and 


da^ 

“dT 


(3.2) 


where the longitudinal field sites z refer to a position in the optical field envelope. 
When the electrons pass through the undulator they exchange information with the 
optical wave envelope in a range of sites. The optical light traveling at speed c, 
remains fixed in z. The slower electrons slip back to site z -x. Those electrons at 
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site z interact with a range of sites in the optical wave envelope. The electrons then 
pass information from one site to another. 

Longitudinal multimode simulations shown in Fig. 3-8 and 3-9 using the same 
parameters as in Fig. 3.7 and 3.8, study the coherence development of the optical 
field. At the bottom-left of each simulation we see the gain over the undulator length. 
The bottom-middle is the single-mode gain spectrum as a function of frequency 
centered around resonance. The bottom right is the evolution of power developed 
through the undulator length, measured in a dimensionless power, corresponding to 
tens of GW. The middle-left plot is the growth in the optical field amplitude over two 
slippage distances in length, -1<z<1. The middle and top-middle show the 
evolution of the power spectrum. The triangle tick mark indicates resonance, v = 0 
and the rectangle tickmark indicates v of the electron beam at the end of the 
undulator. In the klystron undulator, the final power is much higher. The middle-right 
and top-right frames show the development of the electron phase velocity distribution. 
There is a decrease in the average electron phase velocity which corresponds to 
energy given up to the optical field. In the FEL without a klystron, little energy is given 
up to the optical beam, while the system with a klystron gives up significantly more 
energy. The simulation without a klystron and before saturation resulted in an X-ray 
linewidth of AVA. = 0.1%, where AA/A. = Avl2nN and N is the number of periods in the 
undulator. Adding a dispersive section does not change the linewidth. There is no 
degradation to the quality of the optical beam. 

In these systems the standard deviation of phase velocity was oq = 1.04. As the 
standard deviation of the phase velocity inaeases it makes it difficult for the electrons 
to bunch. Tin phase velocity, oq, must be considered when using a klystron FEL. 
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****** FEL Wrap Evolution ****** 


j=932 v^=0 ajj=0 



Figure 3-8. Longitudinal multi-mode simulation for 40 A X ray FEL. 
(Without dispersive section). 

I ****** FEL Wrap Evolution ****** 


j=932 v^=0 ag=0 

0=0.5 N= 413 Og=1.04 Oq=0.26 



Figure 3-9. Longitudinal multi-mode simulation for 40 A X ray FEL. 
(With dispersive section). 
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C. PHASE VELOCITY SPREADS AND KLYSTRONS 


As mentioned in Section II.E, for a change in initial energy Aymc^, there is a 
corresponding change in electron phase velocity, Av = ^nNAyly, in which the Gaussian 
distribution of energies with an rms spread . Ayly, characterizes the standard deviation 
oq = AicA/Ay/y. When oq ~ k, bunching becomes impaired. The use of a klystron to 
bunch the electrons may not have a positive effect if the electron energy spread is so 
large that the modulator cannot develop the phase velocity differences to prepare the 
electrons for bunching. Figures 3-1 Oa, 3-1 Ob, 3-1 Oc, and 3-1 Od are plots on how the 
standard deviation in energy spread can effect the use of a klystron undulator at three 
different dimensionless current. 

Figures 3-1 Oa, 3-1 Ob, and 3-1 Oc show what happens with various energy 
spreads on high-current systems (y = 10000, j = 2200, and y = 500). Figure 3-1 Od is 
the plot for low-current system (y = 1). 

In all four figures you can see as the energy spread, oq, increases, the overall 
gain of the system decreases. By using a klystron undulator it is hoped to increase 
the gain by at least double the magnitude compared to the case where a klystron was 
not used. With a high oq, it is apparent from the figures that a klystron will not 
produce a gain significant enough to warrant its use. The use of a klystron at the high 
oq (2.0 and 3.14) may even cause a decrease in gain. In the monoenergetic beam 
(oq = 0.0), there is no limit to the strength of the dispersive section or to the gain of 
the system; but, this is unrealistic. The plots indicate that a moderate spread, 
0.0 < oq < 1.5, can be compensated for in a klystron and provide a significant gain 
that would warrant its use. The magnitude of D is critical though. If D is to large it 
may have a negative effect and cause the system to decrease in gain. It is also 
shown that as the dimensionless current is increased the critical size of D also 
increases. A lower y can mean a lower D. the distance for the electrons to travel to 
bunch adequately is smaller in the lower energetic beam. 

It was noted that for the 4.5 A and the 1.5 A lasers that SLAG has proposed the 
initial electron phase velocity has a standard deviation of oq = 2.2 and 4.2. The 
strength of the dispersive section thus becomes critical. 
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Figure 3-1 Ob. Gain vs. D plots for dimensionless current j = 2200. 
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D. 1 5 A FEL WITH A KLYSTRON 


The use of a dispersive section to shorten the length of a high power FEL such 
as proposed at SLAG would reduce the cost of the system. The linac is already in 
place, therefore the cost to build the FEL is largely in the manufacturing of an 
undulator. Shortening the undulator would reduce its cost significantly since the 
manufacturing ofundulator sections is large. 

To achieve significant gain in the 1.5 A FEL, the strength of a dispersive section 
for j = 3121 and oq = 4.2 is estimated to be 0.2 from Fig. 3-1 Ob. With high current, 
the maximum gain of an FEL is obtained when vq is slightly off resonance. When 
adding a dispersive section to the undulator, the maximum gain is shifted to a higher 
value of Vq, where vq = nl2D. The optimum for this system is found to be vq ~ 4.0. 

Figure 3-11 is the phase-space plot for a 1.5 A FEL with a length of 30 m, a 
dispersive sfrength of D = 0.2, j = 2170, oq = 3.8, and vq = 4. The length of the FEL 
was shortened from 35 m to 30 m which in turn lowered j and a a ■ The phase-space 
plot shows the system saturating at the end of the undulator. The peak power out is 
14 GW, which is slightly higher than the simulation run without a dispersive section. 
All indications show that a dispersive section can be effective in shortening the length 
of the undulator and increasing gain. 

In order to examine the coherence development, we look at a longitudinal multi- 
mode simulation in Fig. 3-12. We can see in the middle window a narrow spike 
indicating a narrow band width of AAyA,= .1%, which indicates good coherence 
development. The right window shows an average decrease in electron phase velocity 
which corresponds to the electrons giving up energy to the optical beam, creating high 
optical power. 

All the simulations that were run show that a dispersive section may have a 
positive affect in shortening the length of the undulator. The use of a dispersive 
section may decrease the cost of the FEL to be built by SLAG. 
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Figure 3-11. Phase-space plot for 1.5 A system, undulator length of 30 m. 
(With dispersive section) 
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Figure 3-12. Longitudinal multi-mode simulation for 1.5 A X ray FEL. 
(With dispersive section). 
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IV. CONCLUSIONS 


The theory that has been presented is not limited to X ray FELs, but to all FELs. 
The Free Electron Laser has many uses in the military and industrial applications; 
ranging from high-power, speed-of-light self-defense weapon systems to DNA-pair 
imaging. The ability of an FEL to be tunable allows it to overcome limitations that 
other high powered lasers face, such as atmospheric effects. The FEL as a self 
defense weapon system does not require replenishment from a limited magazine. The 
FEL is a system that is needed in the present day military to overcome conditions that 
have outgrown systems now in use. There is a need for the United States to be the 
leader in laser technology to ensure its capability to improve world conditions. 

The proposed SLAC X-ray FEL, which will come on-line in the near future, will be 
the first coherent X ray light source. The acquired technology is not limited to the 
SLAC proposal, but can be used to devebp other types of high powered FELs. An 
injector system now under development for the SLAC linac will enable it to produce a 
4.5 A or 1 .5 A FEL. 

The use of a dispersive section will aid in shortening the length of the undulator. 
A significant part of the cost to build the FEL at SLAC is manufacturing the undulator, 
so shortening of the undulator is economically important. However, beam quality must 
be high in order to use a dispersive section. So the cost of improving the injected 
electron beam quality must be considered in the decision to use a klystron undulator. 
The better the quality of the electron beam the greater the positive affects a dispersive 
section may have on increasing gain and shortening the overall length of the system. 

Finally, we have presented simulation results to show that use of dispersive 
sections is a possible and potential cost-saving option. Further investigation in the use 
of dispersive sections on high-gain FELs is required. For example, positioning the 
dispersive section along the undulator to optimize its effects may be studied. 
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